Reorganization of columnar architecture in the growing visual cortex 
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Many cortical areas increase considerably in size during postnatal development, progressively 
displacing neuronal cell bodies from each other. At present it is unclear to what extent this requires 
remodeling of neuronal circuits. Here, in acute and chronic experiments, we study the layout of 
ocular dominance columns (OD) in cat primary visual cortex (VI) during a period of substantial 
postnatal growth. We find that OD columns reorganize such that the spacing between columns 
is largely unaffected, despite a considerable size increase of VI. We show that the observed mode 
of reorganization can be explained by a fundamental mechanism of growth induced reorganization 
inherent to a very general class of models of activity-dependent formation of OD columns. Our 
work suggests that cortical circuits remain plastic for an extended period in development in order 
to facilitate the modification of neuronal circuits to adjust for cortical growth. 



I. INTRODUCTION 

The brain of most mammalian species grows substan- 
tially during postnatal development. The human brain, 
for instance, weighs on average 350 g in newborns and 
1400 g in adult males [Ij. In cat, the neocortical volume 
increases from ^ 1000 mm^ at birth to ~ 4500 mm^ in 
adulthood [2\. Consistently, the surface area of cat pri- 
mary visual cortex (VI) increases postnatally by 130% 
between week 1 and week 15 [3]. The number of neurons 
remains largely constant during this period [4]. Most of 
the area increase is due to the generation of glial cells, 
the addition of more vasculature and connective tissue, 
and the myelinization of axons; to a lesser extent it also 
reflects the outgrowth and elaboration of axonal and den- 
dritic processes [5J. 

It is likely that such substantial cortical growth in- 
duces remodeling of neuronal circuits. The size increase 
of cat VI implies that the distance between any two neu- 
ronal cell bodies grows on average by ^ 50% during early 
postnatal development. Consequently, synaptic connec- 
tions between cortical neurons that are established early 
in development must be either replaced, or dendritic and 
axonal processes must grow by a similar factor in order 
to preserve these connections. However, such morpho- 
logical changes, albeit subtle, can modify the temporal 
properties of neurons. For instance, signiflcant dendritic 
growth of a pyramidal neuron in cat visual cortex could 
induce delays in dendritic integration [6j on the order 
of a millisecond assuming a lateral dendritic spread of 
^ 300 /im [7J and propagation speeds of voltage changes 
on the order of 100 /j^m/ms [8]. Moreover, since synap- 
tic plasticity can be sensitive on a millisecond time scale 
[9], this might have long lasting eflFects on synaptic ef- 
flcacy. These growth-induced changes potentially apply 
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to a large fraction of synaptic connections and thus may 
aflFect signiflcantly patterns of cortical activity. 

Much of the growth of VI takes place during a period 
when many neurons have already reached fairly mature 
levels of selectivity. For instance, the selective response 
of visual cortical neurons to inputs from one eye or the 
other, called ocular dominance (OD), can already be vi- 
sualized at postnatal week 2 in cat VI [lOj. OD is or- 
ganized into columns which can be labeled within the 
entire VI as early as week 3 [TTJ |12j. Large parts of 
the visual fleld are already represented at this time fT3] . 
Subsequently, between week 3 and week 10, cat VI in- 
creases its size by 50% [12 , 14J, a result which we conflrm 
in Figure [l] In this work, we exploit the fact that OD 
columns can be visualized over a period of signiflcant 
cortical growth to reveal signatures of growth induced 
remodeling of neuronal circuits. 

We analyze the layout of OD columns during corti- 
cal growth in experiment and in models for the activity 
dependent formation of OD columns. Complete recon- 
structions of OD columns in cat VI are obtained using 
2-[^^C]-deoxyglucose (2-DG) and pH]-proline labelling of 
OD columns. Additionally, we follow the development 
of columnar layouts in individual animals by perform- 
ing chronic optical imaging experiments. We flnd that, 
long after OD columns have segregated, columnar lay- 
outs reorganize throughout VI. Despite the considerable 
area increase of VI, the spacing of OD columns is largely 
preserved. To understand these empirical observations, 
we study growth induced reorganization theoretically in 
models for activity-dependent formation of OD columns. 
We argue that in a very large class of models cortical 
growth induces a characteristic zigzag-like bending of OD 
columns by which their spacing is largely restored. This 
class contains, to our best knowledge, any model that 
has been proposed thus far (e.g. [HI [H [HI [H [19]). As 
a representative of this class, we study in detail the re- 
organization of OD layouts in different growth scenarios 
in the Elastic Network model [I9l |20]. We flnd that for 
wide range of growth rates, the growth induced reorga- 
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nization in the model shows the same characteristics as 
the reorganization observed in experiment. Our results 
suggest that cortical growth accompanied by significant 
displacement of neuronal cell bodies can induce substan- 
tial remodeling of neuronal circuits. Furthermore, we 
suggest a possible new role of the critical period in nor- 
mal development which in cat visual cortex overlaps with 
the period of postnatal growth (Fig. : The enhanced 
degree of plasticity during the critical period may facil- 
itate the modifications of neuronal circuits to adjust for 
cortical growth. 



II. RESULTS 

OD Column Spacing Changes Only Little During 
Postnatal Growth of Primary Visual Cortex 

We first measured the size increase of cat primary vi- 
sual cortex (VI) during early postnatal development. We 
labelled complete layouts of OD columns in VI visualized 
by either 2-DG or proline in kittens at different ages be- 
tween postnatal day (PD) 30 and PD84 (N = 18). The 
left panel of Figure [T^ displays a representative 2-DG 
autoradiograph of a kitten at PD30. The right panel dis- 
plays proline labeled OD columns of a kitten at PD60. 
The range of our entire data set compared with the time 
line of cat postnatal development is displayed in Figure 
[T)d. VI was readily discernable by its distinctive colum- 
nar activation pattern in comparison to the labelling in 
surrounding cortical areas (see also [21, 22j). We ob- 
served a size increase of about 30% between two groups 
centered at PD30 and PD70 (r = 0.6248, p < 0.006) 
(Fig. 1|d). To minimize possible infiuences of genetic 
variability [23j, we analyzed a littermate couple on PD30 
and PD72 respectively (green crosses in Fig. [l]b). Con- 
sistently, their cortical sizes also differed by an even larger 
factor of 1.46. Hence, our analysis confirmed the consid- 
erable size increase of VI during early postnatal develop- 
ment observed in previous studies [3| [T2]. 

Given this size increase of VI during the first postnatal 
weeks, we next asked whether the spacing of OD columns 
increases by a corresponding factor over this time period. 
Such a finding would be consistent with a scenario that 
involves no change of column number or layout (some- 
times referred to as "balloon effect" [l2l^4j). First, we 
measured the mean column spacing A of 2-DG/proline 
labehed OD patterns in N = 41 hemispheres (2-DG, N 
= 37; proline, N = 4) between PD28 and PD98 (data 
includes the N = 18 hemispheres used for the analysis 
of VI sizes). To obtain accurate estimations of column 
spacings, we used a wavelet-method that exhibits a pre- 
cision of 2% (relative error, see [23, 25j and Methods). 
As shown in Figure [2^, mean column spacings varied 
between 1.05 mm and 1.28 mm, but did not show a sig- 
nificant increase over this period (r = —0.034, p < 0.825). 
Consistent with this general trend, the column spacings 
of the two littermates differed by less than 10%, despite 



a size difference of almost 50%. 

To follow the development of OD column spacing in in- 
dividual animals we visualized OD columns in cat VI by 
chronic optical imaging experiments over the same period 
(N = 3 hemispheres, total range: PD30-PD94). Figure 
[2Jd shows the layout of OD columns at different ages in 
one of the studied cases (data are high-pass filtered, see 
Methods). Again, we quantified column spacings of OD 
columns by the wavelet-method. As shown in Figure [2]:!, 
we found no systematic increase of column spacings in in- 
dividual animals, thus confirming the conclusions drawn 
from the 2-DG/proline data. Column spacing estimates 
based on optical recordings exhibited a larger variabil- 
ity compared to the 2-DG/proline data. This might be 
explained by the substantial interareal variability of OD 
column spacing [25j and the fact that the recorded re- 
gions spanned only a small portion of VI and possibly 
shifted with age. Taken together, both the 2-DG/proline 
data and the chronic optical recordings demonstrate that 
the postnatal growth of the primary visual cortex is not 
accompanied by a corresponding increase in the spacing 
of OD columns. 



A Universal Mode of Reorganization 

At first glance, the finding of postnatal cortical growth 
without a corresponding increase in OD column spacing 
might seem contradictory. However, we will show in the 
following sections that both findings can naturally be rec- 
onciled within a very large class of models for activity- 
dependent formation of OD columns. 

In most models (e.g. [HI [H [13 [H [19]), the segrega- 
tion of OD columns is caused by a Turing mechanism [26]. 
also called finite wavelength instability [27]. Typically, 
the Turing mechanism results from Mexican-Hat-like cor- 
tical interactions (whether assumed explicitly or arising 
effectively from the assumed dynamics). By this mecha- 
nism, a layer of nonselective neurons segregates into al- 
ternating domains of roughly similar size selective for one 
eye or the other. The spacing between these domains is 
approximately twice the distance between the maximum 
and the minimum of the Mexican-Hat. 

How does cortical growth affect columns that already 
have segregated, according to these models? One possi- 
bility would be that the size of the effective Mexican-Hat 
scales with the size of the cortex and increases by the 
same factor. This could keep neuronal interactions un- 
changed, but would also increase column spacing by the 
same factor. Since we found no evidence for such an in- 
crease of column spacing, we conclude that this scenario 
is rather unlikely. 

Suppose that instead the effective Mexican-Hat does 
not increase or increases less than the cortical tissue. 
Here we conjecture that in this case a reorganization of 
OD domains takes place that shows signatures of the so 
called zigzag (ZZ) instability. The ZZ instability has been 
observed in many inanimate dynamical pattern form- 
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Figure 1: Postnatal growth of cat primary visual cortex (VI). (a) Two representative examples of OD layouts in kittens with 
different age. Left/right panel shows a 2-DG/proline autoradiograph of kittens at PD30/PD60. Yellow lines mark VI borders. 
Scale bar, 10 mm. (b) Area sizes of VI in kittens between PD30 and PD84. The average size increases by ^30% percent 
between the two groups centered around PD30 and PD70. Red line, linear fit (r = 0.6248, p < 0.006). The two green crosses 
mark values from a littermate couple, (c) Growth of cat VI with representation of the visual hemi-field (schematic; HM, 
horizontal meridian; VM, vertical meridian) and time line of cat OD development. OD segregation can be visualized as early 
as PD15 [lOj, OD critical period for monocular deprivation starts at PD21 [TO]. Our dataset includes 2-DG/proline labelled 
hemispheres from kittens between PD28 and PD98 and chronic optical recordings between PD30 and PD96 
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Figure 2: Spacings of OD columns change only little during growth of cat visual cortex, (a) Mean column spacing A of OD 
columns for N = 41 hemispheres of cat VI between PD28 and PD98 (2-DG, N = 37; proline, N = 4). The two green crosses 
mark values from two littermates. Red line, linear fit (r = —0.034, p < 0.825). (b) Chronic recordings of OD columns in an 
individual hemisphere at PD30, PD42 and PD94 (intrinsic signal optical imaging, N = 3, scale bar 1.5 mm). Note that the 
recorded region is only 4-5 column spacings wide and may have shifted during development, (c) Mean column spacings for the 
case in b (blue circles) and for two other cases (red boxes and green crosses). Column spacings from a are shown in light gray 
for comparison. 



ing systems as, for instance, Rayleigh-Benard convec- 
tion, when these systems experience an instantaneous 
isotropic size increase [27]. The result of a ZZ insta- 
bility can be best understood by considering a simple 
pattern of OD columns consisting of alternating stripes 



(Fig. [3^). After an isotropic area increase, the spacing 
between the OD stripes is larger than the spacing set by 
the Mexican-Hat. Subsequently, the structure rearranges 
in a sinusoidal fashion: A zigzag-like bending occurs by 
which the original spacing is largely recovered, locally. 



4 



Figure [3)d illustrates this in Fourier space. The initial 
stripe pattern is described by a wave ~ e*^"^^^^ with spa- 
tial wavenumber /Cmax = 27r/Aniax- Upon area increase, 
its wave vector differs by an amount Sk < from /cmax, 
i.e. the spacing is larger than Amax- Two new Fourier 
modes grow with wave vector (/cmax + Sk)5t -^qyY repre- 
senting the parts of the OD stripes bended to the left and 
right, respectively (Fig. [3^). The diagram in Figure [s]:! 
shows the region of the ZZ instability as a red shaded re- 
gion in the r-5k plane, where r is the control parameter 
of the system expressing the distance from the primary 
(Turing) instability (for a more detailed description, see 
Methods). The diagram also shows that insertion of 
a new stripe or elimination of an existing stripe requires 
strong area increase or decrease, respectively (the regions 
of Eckhaus (EH) instability). The form of this diagram 
is quite general. It is typical of all rotationally invariant 
systems that exhibit a Turing instability and it does not 
reflect the details of the underlying microscopic equations 

Thus, for all models of OD formation that are rotation- 
ally symmetric and based on a Turing mechanism, we 
expect to observe a reorganization of columnar layouts 
upon isotropic area increase that shows signatures of ZZ 
instability. Most models that have been proposed (e.g. 
HH m [HI EHl PJ) fan into this class. In the following, 
we focus on one of its representatives, namely the Elastic 
Network (EN) model [ElEQlEH]. This model has been 
widely used to describe the formation of visual cortical 
selectivity patterns and here serves as a concrete example 
for studying in detail growth induced reorganization. 

Since in a realistic growth scenario the cortex increases 
continuously in size and the initial OD layout is not a 
simple stripe pattern, the growth induced reorganization 
may be more complex than an ideal ZZ instability. How- 
ever, we expect that the reorganization of OD layouts 
still shows signatures of a ZZ instability such as sinu- 
soidal bending or distortion of local pattern elements. 
For simplicity, we refer to all growth induced changes as 
ZZ reorganization in the following. To capture quanti- 
tatively changes in local map layout both in model and 
experiment we calculate a parameter called bandedness a 
[23j that quantifles the tendency of local pattern elements 
to be composed of parallel stripes (see Methods and 
Fig. |1Q| . Large values of a indicate layouts composed 
of regular stripe like parallel domains, whereas small val- 
ues indicate bended, non-parallel stripes or patches. To- 
gether with the column spacing. A, and the bandedness 
a, we use the number of hypercolumns Nhc = A/ h? in 
an area of size A [25j. The quantity A^ is the size of the 
typical unit cell of an OD map, which we refer to as the 
hypercolumn size. Before studying more realistic growth 
scenarios, we consider a near stripe like OD steady state 
solution for the EN model subject to instantaneous area 
increase and reproduce numerically the stability diagram 
sketched in Figure [3]:!. First, however, we study the EN 
model without assuming any growth. 
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Figure 3: What do models predict about OD development 
in a growing visual cortex? Assuming that cortical interac- 
tions do not scale with growth (Fig. [2k, c), most models 
[151 [161 [TT] 18 , 19j predict a zigzag (ZZ) reorganization that 
induces sinusoidal bending of domains (schematic). This is a 
simple mechanism by which the OD stripes can recover their 
column spacing Ai. (a) During isotropic area increase, an OD 
stripe pattern transiently acquires a spacing A2 = Ai + ^A, 
larger than Ai. By the ZZ reorganization, the system locally 
recovers a scale A3, roughly equal to the initial spatial scale 
Ai. (b) In Fourier-space, the initial OD pattern is centered at 
A^maxX. ^A > equals to a shift 5k < 0. Recovery of the ini- 
tial spacing corresponds to the growth of new Fourier modes 

at (/Cmax + Sk)5i ± Qyy with ((/Cmax + Sk)"^ + 

(c) Regions of linear instabilities of the pattern in the {6k, r)- 
plane (redrawn from [29j). The region of the zigzag instabil- 
ity (ZZ, red shaded) is to the left of the Sk = 0-axis, and 
is bounded by a parabola. EH (yellow shaded) denotes the 
regions of the so called Eckhaus instability (see Methods 
and Supplementary Material) where a new OD stripe is 
inserted {Sk < or withdrawn {Sk > For 

— a/t < Sk < — these two regions overlap (striped re- 

gion). Note that the ZZ instability occurs for an arbitrary 
small area increase. Note also that the ZZ instability is the 
generic behavior for 2-dimensional isotropic dynamics with 
linear Mexican-Hat type interactions under the influence of 
isotropic area increase. 



Development of OD Columns in the Elastic Network 
Model Without Growth 

The Elastic Network (EN) model [H [20l El EO] de- 
scribes the development of OD columns by Hebbian plas- 
ticity mechanisms with soft competition, as well as lateral 
intracortical interactions. In this section, we study this 
model under normal conditions (without growth) and 
characterize its Turing mechanism. 
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In the continuous version of the EN model [28], the 
OD map in VI at a given time t is described by a real 
valued field o(x,t), where x represents the position on 
the cortical surface. Negative values of o(x, t) indicate 
a preference for input from the ipsilateral eye, positive 
values a preference for the contralateral eye. This field 
follows a dynamics 

Sto(x, t) = {[so- o(x, t)] A^(x, S, o(-, t)))g + 7? A o(x, t) , 

(1) 

where 

g-(|s.-x|2 + |s.-o(x,^)p)/2a2 
A^(X, S, 0(-, t)) = J^2^g-(|..-y|2 + |.,-o(y,^)P)/2a2 (2) 

is the cortical activity pattern, a controls the receptive 
field size in the stimulus parameter space, (•) denotes the 
average of the ensemble of visual stimuli {S}, 77 mea- 
sures the strength of lateral interactions and A is the 
two-dimensional Laplacian. Selectivities o(x) are modi- 
fied through the cumulative effect of a large number of 
activity events, evoked by the complete stimulus ensem- 
ble. Visual stimuli S = (s^,So) are pointlike and char- 
acterized by a location and an OD value sq, which 
describes whether the activated units are forced to prefer 
the ipsilateral {sq < 0) or the contralateral {sq > 0) eye 
(see Methods). 

By linear stability analysis, we show that OD columns 
segregate by a Turing mechanism. We linearize eq. ([T]) 
around the homogeneous nonselective state o(x) = (see 
Methods). Due to translation and rotation invariance of 
the model, the eigenfunctions of this linearized dynamics 
are Fourier modes ~ e*^^ with eigenvalues [181 EH] 

X{k) = -1 + (1 - e-^'^') - r]^ (3) 

only depending on the absolute value |k| = k (Fig. 
[4^). The amplitude of any small perturbation contain- 
ing a Fourier mode with wave vector k will evolve ~ 
exp(A(/c)t), and therefore spatial frequencies with X{k) < 
are exponentially damped, whereas those with X{k) > 
grow exponentially (shaded region in Fig. |4^). For suf- 
ficiently small cr, X{k) exhibits a positive maximum at 
r{a,r]) = A(A:niax), where /cmax > 0. Thus, the maximum 
will be at some finite wavelength Amax = 27r//Cniax which 
is characteristic for a Turing instability (see also Supple- 
mentary Material). We refer to r as the control pa- 
rameter. Importantly, this defines the timescale r = 1/r 
on which OD columns segregate and the spatial scale 
Amax which will be the dominant scale of the emerging 
pattern of OD columns. Since experiments in cat suggest 
that the segregation of OD columns takes place within a 
few days [10, llj, one may identify r with approximately 
one day in cat postnatal development. Thus, in the EN 
model, OD columns segregate because the nonselective 
state o(x) = becomes unstable and Fourier modes with 
period ~ Amax grow exponentially if a is below a critical 
value. 



Next, we confirmed this by numerical simulations of 
the full model, eq. ([T]). To approximate the stimulus 
ensemble in eq. ([T]) a large, random sample of point- 
like stimuli was drawn at each time step (see Methods 
for details). Different realizations of OD development 
were obtained by presenting different stimulus samples. 
In agreement with previous work fTO", ^20^, a structure 
strongly resembling the OD patterns observed in cat VI 
emerges after a few time steps (see Fig. [l^. Fig. ^p)- 
As a measure of the average cortical selectivity, we mon- 
itor the mean field amplitude A{t) = J (i^x |o(x, t)|. It 
shows a rapid increase during the first lOr and saturates 
thereafter (Fig. At t = 4r the OD pattern already 
exhibits a visible characteristic spacing which becomes 
more dominant in the later time course. Confirming our 
linear stability analysis, we find that the early pattern 
consists of Fourier modes with wavelength ~ Amax (Fig. 
03, d). 

We quantify each simulated OD development by the 
wavelet-method previously applied to our experimental 
data and compute three parameters characterizing the 
map layouts. First, we estimate the mean column spac- 
ing A. We find, that it reaches a constant value after ap- 
proximately lOr close to the predicted value Amax (Fig. 
|4]i). Next, we compute the number of hypercolumns 
Nhc in a simulated area. This number changes only 
very little after OD columns have segregated (Fig. [4^). 
Finally, we calculate the bandedness, a, which quantifies 
the tendency of OD domains to form elongated parallel 
stripes (see Methods for details). In all our simula- 
tions, the initial phase of OD segregation is followed by a 
long stage of slow rearrangements leading to progressively 
more stripe-like OD layouts. This behavior is captured 
by a monotonically increasing bandedness, a, (Fig. |4]f). 
In fact, we numerically find that OD stripes, for which 
the bandedness a is maximal, are stable solutions of the 
EN model. We do not find any other stable solution. 



Zigzag Reorganization in the Elastic Network Model 

Next, we directly tested our hypothesis that the EN 
model shows a ZZ reorganization upon isotropic area 
increase. As a first step, we considered instantaneous 
changes in domain sizes. We initiated simulations with 
a sine wave pattern ~ sin(/cmax^), a state very similar 
to an attractor of the dynamics (Fig. [5]). After a brief 
relaxation period we abruptly changed the size of the 
simulated system (see Methods and Video 2 in Sup- 
plementary Material) and studied its subsequent dy- 
namics. Figure [5] displays snapshots of a typical exam- 
ple of such a simulation (r = 0.15, r] = 0.025), with a 
moderate area increase {Sk/k^s.^ = —0.15). After the 
size reset at t = lOr, stripes start to bend sinusoidally 
(Fig. [5^, upper row). In the power spectrum, the 
growing of additional Fourier modes with wavenumber 
^ {kma^-\-Sk)5t±qyy can be observed (Fig. [5^, lower row 
and compare to Fig. ^jp) . These features are typical for 
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Figure 4: Development of OD column layout in the Elastic Network model without growth, (a) Eigenvalues of the linearized 
dynamics around the nonselective initial state, o(x) = 0. For r(cr, 77) > 0, Fourier modes with positive growth rate form a band 
around /cmax , corresponding to an annulus in 2D. OD columns with spatial scale Amax = 27r//cmax segregate on a timescale 
T = 1/r with r = A(A:max). (b) Development of OD columns starting from initial condition o(x, t = 0) = (77 = 0.025, 
r = 0.2). At t = 4t, the emerging pattern already exhibits characteristics of OD columns. After lOr, columns start to merge 
and progressively reorganize towards a stripe-like layout, (c-f) Time courses of average OD segregation A (c) and of three 
parameters, characterizing OD layouts during development (N = 50 realizations, gray: individual traces, black: example from 
b, red: mean value): (d) mean column spacing A (e) number of hypercolumns Nhc (f) mean bandedness a. Whereas A, 
and Nhc reach mature levels around lOr and exhibit only little variability, the bandedness a increases during the entire time 
course, expressing the fact that solutions slowly converge towards ideal OD stripes. Furthermore, a-time courses display a large 
variability, capturing differences between individual realizations of OD map development. 



the ZZ instability, thus confirming our hypothesis that 
the EN model shows ZZ reorganization upon isotropic 
area expansions. The mean column spacing A exhibits a 
rapid increase at lOr, induced by the instantaneous area 
increase, and subsequently decreases almost to the initial 
value close to Amax (Fig. [5}d). This expresses the fact 
that the ZZ reorganization restores the column spacing 
locally. Consequently, the number of hypercolumns Nhc 
in the OD map increases (Fig. ^p)- The change in layout 
of OD columns due to ZZ reorganization is well captured 
by the drop in the bandedness a (Fig. [5]i). 

Next, we examined the degree of ZZ reorganization 
under various instantaneous size increases. We explicitly 
traced out the regions of stability sketched in Figure ^ 
for the EN model. We analyzed the time scale rzz on 
which ZZ reorganization develops, relative to the time 
scale of column segregation r, by fitting the ZZ modes in 
the power spectrum with Gaussians and measuring their 
peak heights as a function of time. Figure shows 
a typical example of such a fit (r = 0.15, Sk/k^s,^ = 
—0.15, t = llOr). The central peak originates from the 
initial stripe pattern, whereas both of the side peaks are 
the ZZ modes. To explore the full range of expected rear- 
rangements in our model, we systematically vary the con- 
trol parameter, r, and the total increase/decrease, 5L, of 
the linear extent, L, equivalent to varying Sk/k^s,^ (com- 
pare to Fig. [3]:!, see Methods for details). We observe 
three different types of behavior (Fig. ^j^)'- (i) Even 
slight area increases result in an observable sinusoidal 



rearrangement, which can be identified as the ZZ insta- 
bility. Small increases in area result in reorganizations 
on time scales rzz much larger than r. For considerably 
large increases, time scales rzz as low as 2r can be ob- 
served, (ii) For very large instantaneous area increases or 
small control parameters (grey region in Fig. the ini- 
tial pattern decays exponentially after size reset and OD 
segregation restarts with a pattern of characteristic wave- 
length A (pattern plowing) . (iii) OD stripes are found to 
be stable against contraction of the domain {Sk > 0). 
However, no Eckhaus instability (see Supplementary 
Material) is observed, even in case of considerable area 
decrease (Fig. ^p)- 

In conclusion, this analysis shows that for any arbitrar- 
ily small instantaneous area increase {Sk < 0, small), the 
EN model displays ZZ reorganization, evolving on time 
scales larger than the time scale r of OD column segre- 
gation. This effect is robust against a variation of r and 
is, therefore, the generic behavior of the EN model with 
instantaneous area increase. 



Realistic Growth Scenarios 

What specific conclusion can be drawn from these the- 
oretical considerations that can be detected experimen- 
tally? Evidently, an instantaneous area increase may not 
provide the best approximation of cortical growth. More- 
over, a perfect stripe solution may not be a good approx- 
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Figure 5: ZZ reorganization in the EN model after instantaneous isotropic area increase, (a) Snapshots of a simulation starting 
from a near steady state solution of the EN model, a stripe-like OD pattern (upper row) and corresponding power spectra 
(lower row) (77 = 0.025, r = 0.15). At lOr, the linear extent of the system increased abruptly by 18% ((^/c/Zcmax = —0.15). After 
instantaneous area increase, OD domains bend sinusoidally and additional modes appear at (/cmax + Sk)5i zb qyy in the power 
spectra, (b-d) Time courses of mean column spacing A (b), number of hypercolumns Nhc (c), and mean bandedness a (d). 
At lOr, A increases at first by ^ 18%, but subsequently decreases to almost its initial value. Consequently, Nhc increases over 
the considered period, a shows a characteristic drop, reflecting the bending of OD domains. Note that upon area increase the 
EN model displays all features typical for ZZ reorganization. 
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Figure 6: ZZ reorganization in the Elastic Network model under instantaneous isotropic area increase: Timescales and parameter 
regimes. Compare to the general stability diagram sketched in Fig. ^j^. (a) Measuring growth rate of zigzag modes. Black 
crosses mark the power spectrum P{kx = /cmax, ky). Most power is concentrated around the origin ky = 0, stemming from the 
initial plane wave, and at two secondary peaks that correspond to the emerging zigzag modes (compare to Fig. [5]). Red curve 
is a fit of three Gaussians to the power spectrum (Methods), (b) Ratios between the peak height pc of the central Gaussian 
and height pzz of the side peaks at different time points (black crosses). Red line marks linear fit to the logarithm of pzz/Pc 
in the region of exponential growth. The growth rate of zigzag modes is estimated as the slope. In (c), its inverse, the time 
scale rzz, is compared to the intrinsic time scale r of OD segregation for a large set of {6k, r)-values. ZZ reorganization was 
observed for all simulated r > and Sk < 0. For small increases in area size, reorganization grows on time scales rzz much 
larger than the intrinsic time scale r. For considerably large area increases, time scales rzz as low as 2r are observed. In the 
grey region, the original pattern decays and a new OD pattern segregates (pattern plowing). Since simulation time diverges 
for r ^ (Methods), no simulations were carried out in the white region close to r = 0. 
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imation of the OD layout in cat visual cortex. 

Therefore, we next investigated how the ZZ reorgani- 
zation generalizes to more realistic growth scenarios, a 
cortex that is growing continuously beginning before the 
time OD columns segregate. To this end, we performed 
pairs of simulations, one for which we continuously in- 
creased the domain size by various amounts (growing 
system) and one for which we fixed the domain size after 
lOr (non-growing system). To approximate realistic con- 
ditions, we initiated our simulations with the uniform 
and nonselective state o(x) = and linearly increased 
the linear extent of the simulated regions by 20%, 40%, 
60% and 80% between t = Or and t = lOOr (corre- 
sponding to isotropic total area increases of 44%, 96%, 
156% and 224%). Test simulations with sublinear and 
logistic growth did not lead to significantly different re- 
sults. By measuring average column spacing. A, number 
of hypercolumns, Nhc^ and bandedness, a, we quantita- 
tively characterized growth-induced changes and reorga- 
nizations for complex spatial OD layouts under realistic 
growth conditions. 

Figure [7^ shows snapshots of a typical example of 
these simulation pairs (r = 0.16, r] = 0.025, 60% linear 
extension). Shortly after lOr, the non-growing system 
(upper row) tends to converge towards a stripe-like lay- 
out. In contrast, the growing system (lower row) exhibits 
a much more irregular arrangement of columns over the 
whole time course of the simulation. In fact, the growth 
induced reorganization of OD layouts shows aspects of 
the ZZ reorganization. The column spacing of the grow- 
ing system remains close to the value of the non-growing 
system (Fig. [TjD). At the same time, Nhc increases dra- 
matically for the growing system (Fig. ^^)- Comparing 
the time course of the bandedness a for both simula- 
tions (Fig. [rji) shows that growth induced a drop in 
bandedness. Hence, all three parameters show a simi- 
lar time course as for the stereotypic ZZ reorganization 
(compare Fig. [5}D-d). Thus, also for a realistic growth 
scenario, the growth induced reorganization displays the 
fundamental characteristics of a ZZ reorganization. 

To systematically explore how the reorganization de- 
pends on the growth rate and the control parameter, r, 
we quantify the decrease in bandedness for various r and 
total size increases (50 simulations per data point) (Fig. 
[sj. A polynomial fit to the time courses of the banded- 
ness, a, provides us with two characteristic quantities: 
the size of the bandedness drop, Aa, and the duration. 
At, the time interval between the minimum of the drop 
and the preceding maximum (Fig. 8^). Even though 
the simulations exhibit a large variability, bandedness de- 
crease is consistently observed. Figure [sJd demonstrates 
that the percentage of simulations showing a significant 
bandedness drop (with Aa > 0.05, At > 15r) is largely 
independent of the control parameter r. With consider- 
able growth, this percentage is close to 100% as compared 
to < 6% for non-growing systems. For smaller values of 
r the drop is more pronounced. Interestingly, the aver- 
age drop in bandedness depends only weakly on the total 



area increase (Fig. ^^). Thus, even moderate growth 
induces substantial reorganization of OD layouts as com- 
pared to the non-growing system. Perhaps surprisingly, 
the duration At, i.e. the time for which a monotonically 
decreases, is largest for small total area increases. (Fig. 



Reorganization of OD columns During Growth in 
Model and Experiment 

Finally, we looked for signatures of ZZ reorganization 
in experiment and compared it to the model predictions. 
We quantitatively analyzed the 2-DG/proline autoradio- 
graphs by assessing column spacing A, number of hyper- 
columns Nhc^ and bandedness a, and compared the time 
courses of these values with those obtained in our simula- 
tions. As outlined above, our simulations of OD segrega- 
tion together with experimental data [10l[TT] suggest that 
one model time scale, r, corresponds to approximately 
one day in cat postnatal development. In approximate 
accordance with experimental data from Duffy et al. [3], 
we used 60% linear extension between t = and t = lOOr 
corresponding to an area increase of 156%. Similar re- 
sults were obtained when different growth rates are used 
(see Fig. [s]) or when growth was sublinear or logistic 
(see Supplementary Material). 

Figure [9^ displays time courses of the hypercolumn 
number Nhc in 50 simulations obtained for different 
stimulus sets (see Methods and Video 3 in Supple- 
mentary Material). Nhc is linearly increasing with 
simulation time (increase of 47% between t = 30r and 
t = 70r). Figure [oli shows the hypercolumn number 
Nhc for the 2-DG/proline autoradiographs. The num- 
ber of hypercolumns for kittens at later stages in postna- 
tal development is significantly larger, compared to kit- 
tens at early stages (r = 0.7720, p < 0.0002). At PD28, 
VI contains on average 260 ± 40 (N = 5) hypercolumns, 
around PD72 an average of 319 ± 17 (N = 6) (increase of 
23%) hypercolumns. The increase in Nhc in our simula- 
tions stems from the fact that the column spacing A in- 
creases only transiently during continuous isotropic area 
increase (Fig. [9}d and[5]). Interestingly, this behavior, 
characteristic of ZZ reorganization, appears to be present 
also in the data as a closer look at the 2-DG/proline data 
set reveals (Fig. [9^ and Fig. [2^). Column spacings 
are ^5% larger for kittens of intermediate ages between 
PD40 and PD80 and the time course of column spacing in 
model and experiment appear very similar. Furthermore, 
the bandedness a decreases considerably with age, espe- 
cially between P30 and P60. This decrease in bandedness 
is highly significant (r = —0.5754, p < 6 ■ 10~^) indicat- 
ing that OD columns while largely preserving their initial 
spacing reorganize and develop a more isotropic layout 
over time. As shown in Figure |8j such a drop in a is 
characteristic of a ZZ reorganization. Moreover, as for 
the column spacing A and hypercolumn number Nhc^ 
the change in bandedness as well as the time scale on 
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Figure 7: OD development with growth, (a) Snapshots of OD development in the EN model on a domain growing by 156% 
between t = and t = lOOr (upper row) and with the growth stopped at t = lOr (lower row) (r = 0.16, rj = 0.025). (b-d) 
Time courses of mean column spacing A (b), number of hypercolumns Nhc (c), and bandedness a (d) for the simulations in a. 
Despite its large area increase, the growing system displays only a mild and transient increase in A, implying a strong increase 
in Nhc- Note that also a continuous increase of area induces reorganization of OD layouts that shows the characteristic features 
of the zigzag reorganization. 



which this changes evolves are in good agreement to the 
EN model under a realistic growth scenario. We there- 
fore conclude that OD columns reorganize weeks after 
segregation and that models for their activity-dependent 
formation suggest that this reorganization is induced by 
cortical growth. 



III. DISCUSSION 

In cat visual cortex, cortical columns are present long 
before its postnatal growth comes to an end. How does 
cortical growth affect its columnar architecture? Dur- 
ing growth, the cell bodies of most neurons are progres- 
sively displaced from each other. Since this process in- 
evitably requires the modification of neuronal elements, 
it might also alter considerably the patterns of cortical 
activity. Interestingly, the critical period, during which 
circuits are particularly susceptible to changes of activ- 
ity patterns, overlaps with the period of strong postnatal 
growth. It is an appealing hypothesis that cortical cir- 
cuits remain plastic for an extended period in develop- 
ment in order to facilitate the modification of neuronal 
circuits to adjust for cortical growth. Analyzing the or- 
ganization of OD columns in cat primary visual cortex 
during a period of substantial postnatal growth, we pro- 
vide evidence for reorganization of columnar layouts long 
after OD columns have been segregated. Our theoretical 
analysis shows that the observed reorganization is con- 
sistent with a fundamental mechanism of growth induced 
reorganization that is inherent to a very large class of 



models for activity-dependent formation of OD columns. 
By this mechanism, cortical growth induces a charac- 
teristic zigzag-like bending of OD columns which largely 
restores the spacing of columns. Thus, our empirical find- 
ings together with our theoretical analysis suggest that 
cortical growth accompanied by significant displacement 
of neuronal cell bodies can induce substantial remodeling 
of neuronal circuits. 

a. Reorganization of OD Columns during Postnatal 
Growth We visualized the overall layout of OD patterns 
in VI of kittens between PD28 and PD98 by 2DG/proline 
labelling. To minimize genetic variation we included a 
littermate couple into our dataset. In addition, we per- 
formed chronic optical recordings of intrinsic signals to 
follow OD development in individual animals. We ana- 
lyzed the layout of OD columns with high precision using 
a wavelet method [23, 25j. We observed no systematic 
differences in layout between 2DG and proline labelled 
columns [25j. In accordance with [3] |12j [14] we found a 
significant postnatal growth of VI of about 30% over the 
considered period (Fig. [l]), but no corresponding in- 
crease in column spacing A (Fig. [2]). This implies an 
increase in the number of hypercolumns Nhc per visual 
cortex, where a hypercolumn is estimated as a cortical 
region of size A^. Since an attachment of columns at 
the peripheral regions of VI appears unlikely (see Sup- 
plementary Material), these findings suggest that the 
layout of OD columns changes during normal develop- 
ment. In fact, as a signature of ongoing modification of 
columnar architecture during postnatal development, we 
found that the bandedness of OD columns, a, displays a 
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Figure 8: Comprehensive analysis of growth-induced reorga- 
nization of OD column layouts in the EN model, (a) Two pa- 
rameters were used to describe handedness time course (Fig. 
[7|: Size of handedness drop Aa (difference hetween the maxi- 
mal handedness and the minimal suhsequent handedness) and 
duration At (time interval hetween these two handedness ex- 
trema). Blue curve is a polynomial fit to the handedness time 
course (see Methods). Other parameters as in Fig. [t] (b) 
Percentage of simulations with continuous isotropic area in- 
crease at different rates hetween t = and t = lOOr and 
different values of r (N = 50 simulations for each data point, 
error hars indicate s.e.m.) that show a significant handedness 
drop (for which Aa > 0.05 and At > 15r). For considerahle 
growth, this percentage is close to 100% as compared to < 
6% for non-growing systems, (c-d) Average size Aa (c) and 
duration Aa (d) of significant handedness drops. For smaller 
values of r the drop is more pronounced, hut no clear de- 
pendence on the total area increase was found. At is largest 
for relatively small total area increases. Thus, even moderate 
growth induces a substantial drop in handedness a in most 
simulations. 



strong negative correlation with age. Early OD layouts 
exhibit a stronger tendency to form elongated parallel 
bands and hence are qualitatively different from mature 
ones. This decrease in handedness evolves over a period 
of weeks, long after OD columns have segregated. 

b. A Universal Mode of Reorganization Intrigued by 
these findings, we studied as a model for the activity- 
dependent development of cortical circuits the two- 
dimensional Elastic Network (EN) model. Investiga- 
tion of a simple toy growth scenario, OD stripes sub- 
jected to an instantaneous area increase, revealed that 
the generic mechanism of growth-induced reorganization 
in this model is the zigzag reorganization (Fig. [s]). OD 
stripes are found to approximately maintain their intrin- 
sic spacing by bending zigzag-like in response to cortical 
growth (Fig. [s]). This bending occurs on timescales 
Tzz much larger than the time scale r of OD segrega- 
tion. Thus, assuming the latter takes a few days, the ZZ 
reorganization evolves over a few weeks. 



Furthermore, we performed a quantitative analysis of 
OD development in the EN model for realistic growth sce- 
narios. To characterize growth-induced reorganizations, 
the time course of local column spacing A, number of 
hypercolumns Nhc ^ind handedness a were compared to 
non-growing systems. We observed three features charac- 
teristic for zigzag reorganization: (i) a mild and transient 
increase in column spacing A, (ii) a monotonic increase 
in hypercolumn number Nhcj (iii) decrease of local 
anisotropy or handedness a (Fig. [7|[8]). All three fea- 
tures are consistent with our experimental observations 
(Fig. [9]). Indeed, assuming one day in development is 
approximately equal to r, the time scale of OD segre- 
gation in the model (since OD columns are observed to 
segregate within a few days [OT, TT]), the time course of 
growth induced reorganization in the model as captured 
by these parameters exhibits the same characteristic fea- 
tures as the reorganization we observe in experiment. 
Thus, our theoretical analyses indicate that OD columns 
may gradually reorganize by a ZZ mechanism in response 
to cortical growth during postnatal development. 

Importantly, we have argued that ZZ reorganization 
as a mechanism for growth induced reorganization ap- 
plies not only to the EN model, but generalizes to a 
large class of models for the formation of OD columns, 
containing e.g. [IH [H [iTl [H [19]. The ZZ instability 
belongs to the class of so-called "long wavelength insta- 
bilities" [271 ES], describing reorganizations of patterns 
on scales much larger than the spacing of the pattern 
itself (see Fig. ^p)- Thus, although the changes in neu- 
ronal morphologies across VI during postnatal growth 
may be rather subtle on the level of individual neurons, 
they might cause a remodeling of OD columns which is 
coordinated over large parts of VI. 

c. Relation to Previous Work Several experimental 
studies in different mammahan species have shown that 
the cortical surface area as well as the volume increases 
significantly during postnatal growth [3l [121 El 138] . 
Only very few studies have examined the relation be- 
tween postnatal area increase and the spacing of cortical 
columns, and have come to contradictory conclusions. In 
one study [12j, 2-DG labelled OD columns in cat visual 
cortex between 3 and 10 weeks were analyzed. In line 
with our results, this study finds no correlation between 
column spacing A and age, despite a considerable size in- 
crease of VI. Focusing on an earlier period, a longitudi- 
nal optical imaging study of OD columns in single strabis- 
mic cat found an increase of OD column spacing between 
PD27 and PD61 [24J consistent with the slight increase 
we observe over this period (Fig. |9^). The only theo- 
retical study of OD development during cortical growth 
09j focuses on the ID Swindale model [15j. This study 
predicts a splitting of OD stripes in response to cortical 
growth, analogous to the Eckhaus instability discussed 
in the Methods and Supplementary Material. How- 
ever, as acknowledged by the authors themselves, gener- 
alizations to two dimensions cannot be derived from this 
study. Indeed, we show that two-dimensional models for 
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Figure 9: The reorganization of OD layouts in the EN model under a realistic growth scenario is consistent with the reorga- 
nization observed in experiment, (a-c) Model: Time courses of hypercolumn number Nhc (a) column spacing A (b), and 
bandedness a (c) {N = 50 realizations, isotropic total area increase by 156% between t = and t = lOOr, r = 0.16, rj = 0.025). 
Black curve shows the average time course, grey shaded region the SEM. Note the transient increase in column spacing A 
and the typical bandedness drop Aa between 30t and 60t, both characteristic features of the ZZ reorganization, (d-f) Data: 
hypercolumn number Nhc (d), column spacing A (e, reproduced from Fig. |2^), and bandedness a (f) for the same 41 hemi- 
spheres. Nhc increases significantly between the two groups, centered at PD30 and PD70 (r = 0.7720, p < 0.0002). Column 
spacings appear to increase transiently at intermediate ages. The bandedness decreases from an average of 0.15 ±0.04 (N = 8) 
at PD28 to an average of 0.083 i 0.02 after PD84 (N = 6) and exhibits strong negative correlation with age during this period 
(r = —0.5754, p < 6 ■ 10~^). Thus, assuming r ?^lday (see text), the growth induced reorganization in the model as captured 
by these parameters evolves on a similar time scale and exhibits the same characteristic features as in experiment. 



OD development behave qualitatively differently. In par- 
ticular, we found no splitting of stripe-like OD domains 
for isotropic area increases. 

d. Remodeling of Ocular Dominance The class of 
models we propose in this work to explain our empirical 
observations predicts that OD shifts during normal post- 
natal development in a sub-population of neurons. Such 
shifts may require a maximal LGN terminal arbor shift 
of up to A/2 600 /im within one time scale of ZZ reor- 
ganization, corresponding to a speed of shift of ^ 60 /im 
/ day. An important question is whether the observed 
structural changes in in vivo plasticity of geniculocortical 
afferent s are sufficient to account for this reorganization. 
In the normal adult mouse neocortex, structural plastic- 
ity due to elongation and retraction of axonal arbors has 
been observed at a speed of up to 35 fim / day [31]. So 
far, rapid remodeling of LGN afferents in cat visual cor- 
tex in response to short-term monocular occlusion has 
been demonstrated to occur at velocities up to 2 mm / 
day [32j. This suggests that the morphological changes 
required by the mode of reorganization described in this 
paper appear to fall within the range of already observed 
in vivo plastic changes. 

An alternative conceivable scenario might be that the 
observed changes in OD layout during growth is due to a 
particular pattern of displacement of neurons while their 
OD is held fixed. This displacement should then be such 
that the spacing between columns changes only slightly 
and the bandedness decreases progressively. This could 
be achieved, for instance, by a displacement in the form 
of a ZZ-pattern similar to the pattern of reorganization 



of OD in the studied model class. However, it is unclear 
at present what mechanism could give rise to a collective 
displacement of neurons consistent with our empirical ob- 
servations. 

e. A Further Role for the Critical Period According 
to our theoretical analyses, OD development in the Elas- 
tic Network Model can be divided into three phases: (i) 
the spontaneous emergence of spatially isotropic OD pat- 
terns, occurring on a characteristic time scale r {t = — 
lOr) (ii) growth-induced reorganization (t = 10 — lOOr) 
and, (iii) a subsequent phase during which the segregated 
OD patterns slowly reorganize towards stripe-like OD 
layouts {t > lOOr). This third phase covers a period 
in cat postnatal development which is not in the range 
of our data set. To the best of our knowledge, no OD 
layout similar to the stripe-like stationary states of the 
EN model has ever been reported in cat VI. Thus, the 
EN model might not be an appropriate model to describe 
the layout of cat OD columns during this developmental 
stage. The first phase is fairly consistent with experimen- 
tal observations on OD segregation [lOl |33] and suggests 
that r can be approximately identified with one day in 
cat postnatal development. Hence, the second phase of 
growth-induced reorganization in the EN model corre- 
sponds to the phase in cat development during which 
OD columns reorganize. Moreover, this phase largely 
overlaps with the critical period, as depicted in (Fig. 

The impressive ability of cortical circuits to reorga- 
nize during the critical period has been demonstrated 
in numerous studies by artificially manipulating cortical 
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activity, e.g. by monocular deprivation f3i l \35\ \36\ \37]. 
Our findings suggest a possible new role of critical period 
plasticity in normal development. Neuronal cell bodies 
become progressively displaced from each other during 
growth, which inevitably requires morphological changes 
of neuronal elements. Intracortical synaptic connections 
that are established early in development must be either 
replaced or the respective dendritic and axonal elements 
must be elongated. However, this may in turn modify the 
spatial and temporal properties of dendritic integration 
and signal transduction between neurons. Since growth 
induces displacement in a large fraction of pairs of neu- 
rons, this could have a severe effect on cortical processing. 
Critical period plasticity may serve to compensate for 
these growth induced modifications of neuronal elements 
and to facilitate the modifications of neuronal circuits 
necessary to adjust for cortical growth. 



IV. METHODS 
Optical Imaging Data 

Surgery. Anaesthesia was induced with an intra- 
muscular injection of ket amine (10 mg / kg Ketan- 
est(R), Parke-Davis, Berlin, Germany) and xylazine hy- 
drochloride (Rompun(R), Bayer AG, Leverkusen, Ger- 
many) and maintained throughout the experiment using 
nitrous oxide/oxygen anaesthesia (50% N2O/50% O2), 
supplemented with halothane (0.8-1.2%, Eurim Pharma, 
Germany). The EGG, pulmonary pressure, end tidal 
CO2 (3-4%), and rectal temperature (37-38°) were con- 
tinuously monitored. The animals head was fixed in a 
stereotactic frame by means of a metal nut cemented to 
the skull. For optical imaging of VI a craniotomy was 
performed centered at Horsley- Clarke coordinate P4. All 
experiments were performed when the animals were be- 
tween 28 and 94 days old. Successive experiments in 
the same animal were performed with intervals of 7 days 
and in one case 42 days. All animal experiments have 
been performed according to the German Law on the 
Protection of Animals and the corresponding European 
Communities Council Directive of November 24, 1986 
(86/609/EEC). 

Visual Stimulation. Animals were stimulated 
monocularly with high-contrast square-wave gratings 
(subtending 90°x60° visual field) of four orientations (0°, 
45°, 90°, and 135°) moving at a speed of 2cyc/s with a 
spatial frequency of 0.5cyc/deg. Stimuli were generated 
by EZV-Stim software (Optical Imaging Inc., Rehovot, 
Israel) and presented on a LG Electronics Flatron 295 
LCD-monitor (luminosity 180 cd/m^; contrast 300:1; re- 
fresh rate 85 Hz; resolution 1600x1200 pixel) at a dis- 
tance of 25 cm. The eyes were treated with atropine and 
Neosynephrine(R) and refracted appropriately using cor- 
rective corneal contact lenses with artificial pupils with 
a diameter of 3 mm. 

Data Acquisition. The cortical surface was illumi- 



nated by means of two adjustable light guides attached to 
a tungsten-halogen lamp (Spindler & Hoyer, Gottingen, 
Germany) equipped with interference filters for different 
wavelengths. The vascular pattern of the cortex was visu- 
alized at 546 nm =b 10 nm (green), cortical activity maps 
at 707 nm =b 1 nm (red). During data acquisition of in- 
trinsic signals, the camera was focused 650-750 /im below 
the cortical surface. A tandem-lense was used for imag- 
ing [44J. The ORA 2001 system (Optical Imaging Inc.), 
equipped with a cooled Theta CCD system (384x288 
pixel chip from Thomson-CSF) was used for collecting 
the intrinsic signals. We acquired a series of frames ev- 
ery 12 s, whereby a grating of a given orientation was 
presented for 2s in a static mode, followed by 4.2 s of 
data acquisition during which the grating was moved in 
both directions along the axis orthogonal to its orienta- 
tion. We used episodic stimulation during data acquisi- 
tion (7 frames of 600 ms duration). The first frames were 
excluded from further analysis. The stimulus presenta- 
tion was monocular, an eye shutter was used to conceal 
the eyes. A single stimulus trial consisted of 2x8 stimu- 
lus conditions (4 grating orientations for the left and the 
right eye) and 8 isoluminant blanks presented in a ran- 
dom sequence. Twelve trials were usually presented to 
obtain a map, so that every stimulus was shown 24 times. 
We first calculated "single condition maps" in which the 
images acquired during presentation of a particular stim- 
ulus were divided by the sum of all different stimulus 
conditions ("cocktail blank procedure") [IF, "46]. Differ- 
ential maps for ocular dominance (OD) were calculated 
by summing all left eye activity maps and subtracting all 
right eye activity maps. 

Data Preprocessing. All computed differential maps 
were preprocessed in order to remove overall variations in 
signal strength and measurement noise. To achieve this, 
first the local average signal strength 

was calculated using the Fermi filter kernel K{y) = 
1/ (1 + exp ((y^ — hp)/P)) and subtracted from the pat- 
tern 

/i(x) = /o(x)-/(x) 

We used (3 = 0.2 and a cutoff wavelengths hp = 1.8 mm. 
Then we lowpass filtered the image with a second Fermi 
filter using a cutoff wavelengths of Ip = 0.7 mm and 
(3 = 0.2. The pattern was then centered to yield 
J^^r^ (Px I{x.) = and its variance was normalized to 
one. This overall bandpass filtering ensured that struc- 
tures on a scale between 0.7 mm and 1.8 mm were only 
weakly attenuated by the preprocessing and enabled us 
to do further quantitative analysis. 
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2 DG - Autoradiography 

Animals. We analysed patterns of OD columns in 
primary visual cortex (VI) of 31 cats (49 hemispheres). 
Some of the OD maps of these animals have already been 
published (see [TTJ [25] and references therein). OD pat- 
terns were labelled with either 2-[^^C]-deoxyglucose (2- 
DG) autoradiography \47\ after monocular stimulation of 
the animals or by [^H]-proline autoradiography after in- 
jection of the labelled proline into one eye which labels 
the thalamocortical afferents of that eye in cortical layer 
IV [IHlllS]. All experimental procedures are described in 
detail in [211 [22l \50]. A list of the dataset as well as a 
description of the preprocessing is given in [25j. 



Data Analysis Method 

Column spacing A and handedness a of both, data 
and simulations, were analyzed using the wavelet-method 
introduced in [23j. 

Estimation of Column Spacing 

For each OD pattern /(x) we calculated a wavelet rep- 
resentation, using wavelets which covered only a few hy- 
percolumns but exhibited a strong periodicity. These 
representations were obtained from 
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where x, ^, / are the position, orientation and scale of the 
wavelet ipx,9,i and I{^^0^ I) denotes the array of wavelet 
coefficients. We used complex Morlet wavelets defined by 
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is the two-dimensional rotation matrix. 

The characteristic wavelength of a wavelet with scale 
/ is A^/ with A^ = 27r/|k^|. a denotes the anisotropy 
of the wavelet. We used relatively large (k^ = (7,0)) 
isotropic wavelets to estimate local column spacing A(x). 
We first calculated 
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of the wavelet coefficients for every position x and 
then determined the scale by computing [(x) = 
argmaxi{I{'K^l)). The corresponding characteristic lo- 
cal wavelength was then A(x) = /(x)A^. We used 12 
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Figure 10: Wavelet analysis of columnar layouts, (a) Exam- 
ples of two wavelets superimposed on a stripe-like region of a 
simulated OD pattern. The real parts of the complex-valued 
wavelet are depicted. Positive regions are delineated by red 
lines, negative regions are delineated by green lines. The 
wavelet of optimal spatial frequency and orientation (solid 
line) and one of non optimal orientation (dashed line) are 
shown, (b) Wavelets superimposed on sinusoidally bended 
stripe patterns, characteristic for the ZZ reorganization, (c- 
d) The normalized squared modulus of the wavelet coefficients 
as a function of orientation for the pattern in a (c) and for 
the pattern in b (d). Note that for the stripe-like pattern 
|/(x, ^)|^ is strongly modulated and exhibits a pronounced 
peak. The sinusoidally bended stripes lead to a broader 0- 
dependence and a less pronounced peak. This difference is 
captured by the handedness parameter a. 



equally spaced orientations and 16 scales Ij between 0.5 
mm and 2.0 mm. [(x) was then estimated as the maxi- 
mum of a polynomial of 6th degree in fitting the /(x, Ij) 
(least square fit) for a given position x. Based on the lo- 
cal column spacing A(x) we calculated the mean column 
spacing A = (A(x))^ and the number of hypercolumns 
in a given area A, Nhc = Aj IS? . Here and in the fol- 
lowing, (•)^ denotes the average over VI or the simulated 
domain. 

Estimation of Local Bandedness a The orienta- 
tion dependence of the wavelet coefficients was used to 
calculate a parameter measuring the anisotropy of local 
pattern elements as shown in Fig. |1Q| For a pattern con- 
sisting of elongated bands, the magnitude of the wavelet 
coefficients depends strongly on the wavelet orientation 
and is largest if the orientation of the wavelet matches the 
orientation of the bands (Fig. 10 1, c). For a pattern 



consisting of more bended stripes or isotropic patches, 
the wavelet coefficients depend only weakly on the orien- 
tation of the wavelet (Fig. [lOf p, d). Using only wavelets 
with local wavelength A(x), we therefore calculated 



.'(x) 



£ d0\i{x,0)\^e^^'' ^ £ de\i{x,0)\^ . (4) 
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We defined the local handedness as 



s(x) = / , 

Jvi 



d^yK{^-y)s\y) 



Ivi 



d'yK{^-y) , (5) 



where i^(x) = exp (-x2/2cr2) and a = 1.3A. To 
estimate Q and ^fj), wavelet coefficients were computed 
for 9 equally spaced orientations. With this choice of 
cr, 5(x) is sensitive to the occurrence of band-like regions 
that extend over the size of a hypercolumn. Based on the 
local handedness s(x) we computed the handedness a = 
(^(x))^ to characterize the overall layout of simulated 
and measured OD patterns. 



Parameter Range of ZZ instability 

The stability diagram shown in Fig. ^ is very gen- 
eral, applying to the class of two dimensional translation 
and rotation symmetric systems in which the pattern 
forming process is based on a Turing-mechanism [27 \. 
Assuming a plane wave solution t)e*^"^^^^, close 

to the Turing-instability its behavior is governed by the 
Newell- Whitehead equation [51] 



dtA = rA^ (5^ - idy2fA - \AfA , 



(6) 



were r is called the control parameter of the system. 
The form of eq. (|6| is independent of the microscopic 
details of the instability mechanism [27]. Inserting the 
ansatz A = |^|e<^fca^+0o)^ yields a uniform solution 
\A\ = \Jt — 5k?' as long as \5k\ < ^r. The wavenum- 
ber of this solution differs by an amount 5k from /Cmax, 
i.e., their characteristic spacing is either smaller or larger 
than A. Linear stability analysis yields the behavior of 
this solution depending on the control parameter r and 
the difference in wave vector 5k [52]. Considering general 
perturbations of the form 5 A ^ ^'^QyV^'^Q^x^ ^^le uniform 
solution is found to be unstable against perturbations 
with = as soon as 

5k <0. 

This defines the domain of the ZZ instability (red shaded 
regions Fig. ^p)- In contrast, for Qy = the uniform so- 
lution is only unstable if\5k\ > ^Jr jZ. This is the domain 
of the so-called Eckhaus instability (EH; yellow shaded 
regions in Fig. [3]:!; see Supplementa ry M aterial for 
an illustration). Yox < 5k < -\frjz EH and ZZ 

overlap (striped region in Fig. and model specific 
properties determine which of the two behaviors is ob- 
served in this regime. For \5k\ > y^, eq. (|6| has no 
uniform stationary solution. In this regime, the initial 
stripe pattern decays exponentially and OD segregation 
restarts {pattern plowing). 



Linearization of the Elastic Network Model 

First, we demonstrate that the homogeneous nonselec- 
tive state o^orn(x) = is a stationary solution of eq. ([T]) 



]. Inserting oiiom 

(x) = into eq. ^ yields 



Aa-{'^^ S, Ohon 



. |sr-x| 



27rcr2 



After averaging over the ensemble of stimuli with nor- 
malized uniform densities pr and ps^, we obtain 



0. 



This shows that Ohom{'^) = is a fixed point of the dy- 
namics in eq. ([T]). The stability of this fixed point can 
be studied by linearizing the r.h.s of eq. ([T]) [HI [28], 



J^[o,X,cr, S] ^ J^[Ohom,X,cr, S] + J 



5T 



o(y). 



Assuming a uniform stimulus density ps^ across the cor- 
tical surface this results in 



dto{x.) = [t] A 



1 o(x)- 



I 



dye" 



_ (x-y)f 



o(y)- 



(7) 

Note that the linearization is only governed by the vari- 
ance of the stimulus ensemble, and not by higher order 
statistical moments and is invariant under spatial rota- 
tions and translations. Due to the latter, its eigenfunc- 
tions are Fourier modes ~ e*^^ and their eigenvalues \{k) 
only depend on the absolute value /c = |k|. By Fourier- 
transforming eq. ([7| with (5^) = 1 we obtain \{k) in eq. 

m 

The homogeneous nonselective state Ohom{^) = is 
unstable, if some eigenvalue with /c > is larger than 
zero. For r] > and a > 0, X{k) has a single max- 
imum at /cmax = ^\/ln(l/7^). A(/Cniax) is positive if 
the width a of the activation function is smaller than 



cr* = ^/l — T] -\- r]\nr] . In this regime, the homogeneous 
nonselective state is unstable and the maximum posi- 
tive eigenvalue r = X{kmax) defines a timescale r = 1/r 
on which OD columns segregate and a spatial scale 
Amax = 27r//cniax that approximately corresponds to the 
spacing of columns. 



Numerical Integration 

Simulations were performed on a 64x64 grid with pe- 
riodic boundary conditions. Simulated systems were spa- 
tially discretized with at least 4 grid points per Amax to 
achieve sufficient resolution in space. Test simulations 
with larger grid sizes (128x128, 256x256) did not lead 
to significantly different results. Progression of time was 
measured in units of the intrinsic timescale r. Here, the 
integration time step 5t is bounded by the relevant decay 
time constant of the Laplacian in eq. ([T]) around /c^ax 
and by the intrinsic timescale r of the system. We used 
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6t = min {1/(20/^ /c^g^^), r/lO} to ensure good approxi- 
mation to the temporally continuous changes of the OD 
patterns. Note that, since r = 1/r, simulation time di- 
verges for r ^ 0. The EN model dynamics was simulated 
using an Adams-Bashforth scheme for the first term on 
the r.h.s. of eq. ([T]). The second term was treated by 
spectral integration, exhibiting unconditional numerical 
stability. The stimulus parameters Sr and Sq were cho- 
sen to be uniformly distributed with densities pr and p. 
such that (5^) = 1. The stimulus average in eq. M 
was approximated by choosing a random representative 
sample of Ng stimuli at each integration time step, with 

where n is the dimensions of the feature space (in our 
case, n = 3), = (L/Amax)^ the size of the simulated 
system in units of A^^^^, Sg the resolution in feature space, 
A^o the number of stimuli we required to sufficiently ap- 
proximate the cumulative effect of the ensemble of stim- 
uli within each feature space voxel, and ^00 = ^^^oJe 
a factor that depends on the specific form of nonlin- 
ear competition between Fourier modes in the EN model 
(unpublished data and personal communication from M. 
Huang and F. Wolf), y/rgoo is proportional to the in- 
verse of the expected mean amplitude of the OD pattern 
(unpublished data). With A'o = 50 and Ss = 0.15 we 
ensured a low amplitude to noise ratio for all the simu- 
lations. Typical values for Ng were between 40000 and 
200000. To model development prior to OD segregation, 
we initialized simulations with o(x, t = 0) = 0. Different 
realizations of OD development were obtained by using 
different stimulus samples. 

Instantaneous Area Increases. We started sim- 
ulations with a predefined system length of Lq and 
used a parameter dependent sine-wave initial condition 
sin(A:max x) with periodicity A^ax = 27r/A;max 
and amplitude close the amplitude of the expected sta- 
tionary state. The system was integrated keeping the 
length of the system fixed for the first lOr. Within this 
period the amplitude A{t) = / (Px\o{-K^t)\ of the OD 
stripe patterns relaxed towards a stationary value while 
OD layout underwent only minor random fiuctuations. 
Next we rescaled the system length to Lo-\-SL (no change 
in number of grid points), where 

max^^ 
{Sk/kmsc^) - 1 

is determined from the desired value of Sk/k^s,^, thereby 
increasing {SL > 0) or decreasing {SL < 0) the total area 



size covered by the OD pattern. Accordingly, we ad- 
justed the number of stimuli, Nsj and, since A ~ l/L'^j 
the matrix for the spectral integration step. Note that 
the numerical value of a, i.e. the width of the activa- 
tion function Acr(x, S, o(-, t)), remains unchanged by this 
rescaling, but its extension relative to L decreases (and 
the number of grid points per activation blob decreases 
as well). After rescaling, we continued simulations for 
another 190r. 

To estimate time scales rzz of the ZZ reorganization, 
we monitored the growth of the two ZZ modes over time 
(see Fig. |6^, b). At each time point, we fitted Gaussians 
of variable size, width and position to the original mode 
and the two growing ZZ modes along the expected axis of 
the ZZ modes in Fourier space (least square fit) (see Fig. 
|6^). The timescale of ZZ reorganization was extracted by 
linear fitting (least square fit) the logarithm of the ratio 
between the peak height of the ZZ modes, pzzj and the 
central peak height stemming from the original pattern, 
Pc (see Fig. |6}d) . We constraint the fit to values between 
15% and 60%^of the maximal peak ratio to capture best 
the regime of exponential growth of the ZZ modes. The 
fitting procedure was only applied to simulations with 
maximal peak ratio of at least 0.05. 

Continuous Area Increases. To approximate real- 
istic growth conditions, we linearly increased the linear 
extent L of the simulated regions between t = Or and 
t = lOOr (see also Supplementary Material). System 
length L, the Laplacian A, and the number of stimuli 
Ns were updated at each integration step. After lOOr, 
integration was continued for another lOOr with fixed lin- 
ear extent. To quantify the size of the bandedness drop 
Aa due to growth as well as the time At over which a 
decreased, we fitted an 8th-order polynomial to each a- 
time series (least square fit). This captured the coarse 
grained development but smoothed out fiuctuations on 
intervals smaller than 3r (see Fig. |8^). We only consid- 
ered bandedness decreases of minimum size Aa > 0.05 
and duration At > 15r in our analysis. 

Statistics, r-values denote Pearson's linear correla- 
tion coefficient; p- Values were obtained with Student's 
t-tests. All error bars and ± are s.e.m.. 



Acknowledgements 

We express our gratitude to F. Wolf for inspiring dis- 
cussions. We thank M. Huang and F. Wolf for comments 
on numerical procedures. We thank P. Mehta and S. 
Palmer for helpful comments on the manuscript. 



[1] Dekaban, A.S. & Sadowksy, D. Changes in brain weight 
during the span of a human life: relation of brain weight 
to body heights and body weights. Ann. Neurol 4, 345- 



356 (1978). 

[2] Villabanca, J.R., Schmanke, T.D., Crutcher, H.A. , Sung, 
A.C. , & Tavabi, K. The growth of the feline brain from 



16 



fetal into adult life i. a morphometric study of neocortex 

and white matter. Dev. Brain Res. 122 , 11-20 (2000). [24 
[3] Duffy, K.R., Murphy, K.M., k Jones, D.G. Analysis of 

the postnatal growth of visual cortex. Vis. Neurosci. 15, 

831-839 (1998). [25 
[4] Cragg, B.G. The development of synapses in the visual 

system of the cat. J. Comp. Neurol. 160, 147-166 (1975). 
[5] Purves, D. Neural activity and the growth of the brain. 

(ed. Brozolo, L.) (Cambridge University Press, 1994). 
[6] London, M. & Hausser, M. Dendritic computation. Ann. [26 

Rev. Neurosci. 28, 503-32 (2005). 
[7] Martin, K.A. k Whitteridge, D. Form, function and in- [27 

tracortical projections of spiny neurones in the striate 

visual cortex of the cat. J. Physiol. 353, 463-504 (1984). [28 
[8] Hausser, M., Spruston, N. & Stuart, G.J. Diversity and 

dynamics of dendritic signaling. Science 290, 739-44 [29 

(2000) . 

[9] Dan, Y. & Poo, M.M. Spike timing-dependent plasticity 

of neural circuits. Neuron 44, 23-30 (2004). [30 
[101 Crair, M.C., Horton, J.C., Antonini, A., k Stryker, M.P. 
Emergence of ocular dominance columns in cat visual 
cortex by 2 weeks of age. J. Comp. Neurol. 430, 235-249 

(2001) . [31 
[11] Rathjen, S. k Lowel, S. Early postnatal development 

of functional ocular dominance columns in cat primary 
visual cortex. Neuroreport 11, 2363-7 (2000). [32 
[121 Rathjen, S. Schmidt, K.E., k Lowel, S. Postnatal growth 
and column spacing in cat primary visual cortex. Exp. 
Brain Res. 149, 151-158 (2003). [33 
[131 Sireteanu, R. k Maurer, D. The development of the kit- 
ten's visual field. Vis. Res. 22, 1105-11 (1982). [34 
[141 Rathjen, S. et al. The growth of cat cerebral cortex in 
postnatal life: a magnetic resonance imaging study. Eur. 
J. Neurosci. 18, 1797-1806 (2003). [35 
[151 Swindale, N.V. A model for the formation of ocular dom- 
inance stripes. Proc. Royal Soc. London, B 208, 243-264 
(1980). [36 
[161 Miller, K.D., Keller, J.B. k Stryker, M.P. Ocular dom- 
inance column development: Analysis and simulation. 
Science 245, 605-615 (1989). [37 
[171 Obermayer, K., Blasdel, G., k Schulten, K. Statistical- 
mechanical analysis of self-organization and pattern for- [38 
mation during the development of visual maps. Phys. 
Rev. A 45, 7568-7589 (1992). 
[181 Scherf, O., Pawelzik, K., Wolf, P., k Geisel, T. Theory of 

ocular dominance pattern formation. Phys. Rev. E 59, [39 
6977-6993 (1999). 
[191 Goodhill, G.J. k Cimponeriu, A. Analysis of the elastic 

net model applied to the formation of ocular dominance [40 
and orientation columns. Network (Bristol, England) 11, 
153-168 (2000). 

[201 Durbin, R. k Mitchison, G. A dimension reduction [41 
framework for understanding cortical maps. Nature 343, 
644-647 (1990). 
[211 Lowel, S. k Singer, W. Monocularly induced 2- 

deoxyglucose patterns in the visual cortex and lateral [42 
geniculate nucleus of the cat. IL awake animals k stra- 
bismic animals. Eur. J. Neurosci. 5, 857-869 (1993). 
[221 Lowel, S. k Singer, W. Monocularly induced 2- [43 
deoxyglucose patterns in the visual cortex and lateral 
geniculate nucleus of the cat. L anaesthetized and paral- 
ysed animals. Eur. J. Neurosci. 5, 846-856 (1993). [44 
[231 Kaschube, M., Wolf, P., Geisel, T., k Lowel, S. Genetic 
influence on quantitative features of neocortical architec- 



ture. J. Neurosci. 22, 7206-7217 (2002). 
Miiller, T. et al. An analysis of orientation and ocular 
dominance patterns in the visual cortex of cats and fer- 
rets. Neur. Comp. 12, 2573-2595 (2000). 
Kaschube, M. et al. The pattern of ocular dominance 
columns in cat primary visual cortex: intra- and in- 
terindividual variability of column spacing and its de- 
pendence on genetic background. Eur. J. Neurosci. 18, 
3251-3266 (2003). 

Turing, A. The chemical basis of morphogenesis. Phil. 
Trans. Royal Soc. London, B 237, 37-72 (1952). 
Cross, M.C k Hohenberg, P.C. Pattern formation out- 
side of equilibrium. Rev. Mod. Phys. 65, 854-1086 (1993). 
Wolf, P. k Geisel, T. Spontaneous pinwheel annihilation 
during visual development. Nature 395, 73-78 (1998). 
Manneville, P. Dissipative Structures And Weak Turbu- 
lence, (ed. Araki H., Libchaber A. k Parisi, G.) (Aca- 
demic Press, 1990). 

Erwin, E., Obermayer, K., k Schulten, K. Models of 
orientation and ocular dominance columns in the visual 
cortex: a critical comparison. Neur. Comp. 7, 425-468 
(1995). 

De Paola, V. et al. Cell type-specific structural plasticity 
of axonal branches and boutons in the adult neocortex. 
Neuron 49, 861-875 (2006). 

Antonini, A. k Stryker, M. P. Rapid remodeling of ax- 
onal arbors in the visual cortex. Science 260, 1819-1821 
(1993). 

Crowley, J. k Katz, L.C. Early development of ocular 
dominance columns. Science (2000). 
Hubel, D.H. k Wiesel, T.N. The period of susceptibility 
to the physiological effects of unilateral eye closure in 
kittens.sk J. Physiol. 206, 419-436 (1970). 
Shatz, C.J. k Stryker, M.P. Ocular dominance in layer 
iv of the cat's visual cortex and the effects of monocular 
deprivation. J. Physiol. 281, 267-83 (1978). 
Antonini, A. k Stryker, M.P. Plasticity of geniculocor- 
tical afferents following brief or prolonged monocular oc- 
clusion in the cat. J. Comp. Neurol. 369, 64-82 (1996). 
Hensch, T.K. Critical period plasticity in local cortical 
circuits. Nat. Rev. Neurosci. 6, 877-888 (2005). 
Drenhaus, U., Rager, G., Eggli, P. k Kretz, R. On the 
postnatal development of the striate cortex (vl) in the 
tree shrew (tupaia belangeri). Eur. J. Neurosci. 24, 479- 
490 (2006). 

Oster, A.M. & Bressloff, P.C. A developmental model of 
ocular dominance column formation on a growing cortex. 
Bull. Math. Biol. 68, 73-98 (2006). 

Marotte, L.R. Goldfish retinotectal system: continuing 
development and synaptogenesis. J. Comp. Neurol. 193, 
319-334 (1980). 

Raymond, P. A. k Easter, S.S. Postembryonic growth 
of the optic tectum in goldfish, i. location of germinal 
cells and numbers of neurons produced. J. Neurosci. 3, 
1077-1091 (1983). 

Easter, S.S. k Stuermer, C.A. An evaluation of the hy- 
pothesis of shifting terminals in goldfish optic tectum. J. 
Neurosci. 4, 1052-1063 (1984). 

Constantine-Paton, M. k Law, M.L Eye-specific termi- 
nation bands in tecta of three-eyed frogs. Science 202, 
639-41 (1978). 

Ratzlaff, E.H. k Grinvald, A. A tandem-lens epifluo- 
rescence macroscope: Hundred-fold brightness advantage 
for wide field imaging. J. Neurosci. 36, 127-137 (1991). 



17 



[45] Bonhoeffer, T. & Grinvald, A. Optical imaging based on 
intrinsic signals: The methodology, in Brain Mapping. 
The Methods (ed. Toga, A., Mazziotta, J.) (Academic 
Press, 1996). 

[46] Engelmann, R., Crook, J.M. & Lowel, S. Optical imaging 
of orientation and ocular dominance maps in area 17 of 
cats with convergent strabismus. Vis. Neurosci. 19, 39- 
49 (2002). 

[47] Sokoloff, L. et al. The [14c]deoxyglucose method for the 
measurement of local cerebral glucose utilization: theory, 
procedure, and normal values in the conscious and anes- 
thetized albino rat. J. Neurochem. 28, 897-916 (1997). 

[48] Grafstein, B. Transneuronal transfer of radioactivity in 
the central nervous system. Science 172, 177-179 (1971). 

[49] Wiesel, T.N., Hubel, D.H., & Lam, D.M. Autoradio- 
graphic demonstration of ocular-dominance columns in 
the monkey striate cortex by means of transneuronal 
transport. Brain Res. 79, 273-279 (1974). 

[50] Lowel, S. Ocular dominance column development: Stra- 
bismus changes the spacing of adjacent columns in cat 
visual cortex. J. Neurosci. 14, 7451-7468 (1994). 

[51] Newell, A.C. & Whitehead, J. A. Finite bandwidth, finite 
amplitude convection. J. Fluid Mech. 38, 279-303 (1969). 

[52] Busse, F. & Whitehead, J. A. Instabilities of convection 
rolls in a high prandtl number fluid. J. Fluid Mech. 47, 
305-320 (1971). 

[53] LeVay, S., Stryker, M.P. & Shatz, C.J. Ocular dominance 
columns and their development in layer iv of the cat's 
visual cortex: a quantitative study. J. Comp. Neurol. 
179(1), 223-44 (1978). 

[54] Luskin, M.B. & Shatz, C. Studies of the earliest gen- 
erated cells of the cat's visual cortex: Cogeneration of 
subplate and marginal zones. J. Neurosci. 5(4), 1062- 
1075 (1985). 

[55] Luskin, M.B. & Shatz, C.J. Neurogenesis of the cat's 
primary visual cortex. J. Comp. Neurol. 242(4), 611-31 
(1985). 

[56] Crampin, E.J., Gaffney, E.A., & Maini, P.K. Reaction 
and diffusion on growing domains: Scenarios for robust 
pattern formation. Bull. Math. Biol. 61(6), 1093-1120 
(1999). 



18 



V. SUPPLEMENTARY MATERIAL 
Adding OD columns at the periphery of VI 

At a first glance, a different scenario of postnatal re- 
organization during VI growth seems conceivable which 
assumes the addition of new functional modules at the 
periphery, thereby preserving the early OD layout in the 
central regions. Here we illustrate that this scenario in 
fact implies substantial large-scale rearrangements of OD 
columns and is therefore rather implausible. 

During the period of postnatal growth considered in 
this paper, the number of neurons remains largely con- 
stant [4J. The increase in total area is largely due to 
the generation of glial cells, the addition of more vas- 
culature and connective tissue, and the myelinization of 
axons, but also involves an increase in neuropil, i.e. the 
outgrowth and elaboration of axonal and dendritic pro- 
cesses and extensive remodeling of synaptic connections 
[H m [53l [SH |55] . Figure 11 schematically demonstrates 
that under these assumptions, the addition of new OD 
columns at peripheral regions of VI while preserving the 
OD layout in central regions of VI implies substantial 
remodelling of OD in large parts of the neuronal popula- 
tion. 

Ocular dominance columns can already be visualized 
at postnatal week 2 in cat VI [lOj. This layout of alter- 
nating domains is sketched in the upper row of Figure 
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The lower row displays the potential layout of OD 
domains in the adult animal. Additional columns are in- 
serted at peripheral regions. The column spacing A is 
the same for both OD patterns. Furthermore, the cen- 
tral part of the early and mature OD layout are identi- 
cal (dashed lines). However, neurons labelled with a red 
number in the mature animal have changed their pre- 
ferred eye. Note their wide distribution across the cor- 
tical surface. Note further that during the continuous 
expansion of VI, eye preference changed multiple times 
in individual termination zones, in particular at the pe- 
ripheral regions (e.g. regions 1,18). Thus, we conclude 
that this mode of postnatal reorganization is implausible 
to occur during normal development. 



Eckhaus Instability 

The Eckhaus instability is a mode of reorganization 
predicted to occur for models with effective Mexican-Hat 
interactions in isotropically growing as well as shrinking 
systems. Rearrangements of domains consist either in 
inserting new domains (area increase) or removing ex- 
isting ones (area decrease). Even if it takes place in a 
two-dimensional system, it is effectively one-dimensional 
reorganization. It does not require rotation symmetry of 
the Turing Mechanism. 

The result of an Eckhaus instability can be best un- 
derstood by considering a simple pattern of OD columns 
consisting of alternating stripes (Figure 12). After 



an isotropic area increase, the spacing between the OD 
stripes is larger than the spacing set by the Mexican- 
Hat. Subsequently, new stripes of OD domains emerge 
between existing ones by which the original spacing is 
eventually recovered (Fig. |12^ ). Figure |12[ 3 illus- 
trates this in Fourier space. The initial stripe pattern is 
described by a wave ~ e^^r^^^^ with spatial wavenumber 
^max = 27r/Amax- Upou area increase, its wave vector 
differs by an amount < from /Cmax- Two Fourier 
modes grow with wave vector ±/cmaxX, representing the 
new OD pattern in Figure \12\ » with similar spacing as 
the original pattern. 

Similarly, after an isotropic area decrease, the spac- 
ing between the OD stripes is smaller than the spacing 
set by the Mexican-Hat. Subsequently, a fraction of the 
OD stripes is removed and the original spacing is even- 
tually recovered (Fig. 12:^). Figure [T2] d illustrates this 
in Fourier space. Upon area decrease, the OD pattern's 
wave vector differs by an amount < from /Cmax- Two 
Fourier modes grow with wave vector ±/cmaxX, represent- 
ing the new OD pattern in Figure [12] :; with similar spac- 
ing as the original pattern. As shown in our manuscript 
the domain of the Eckhaus instability in the (r, (5A:)-plane 
is \5k\ > \JVJ^ (see Figure 3c in the manuscript). For 
< 0, this domain largely overlaps with the domain of 
the ZZ instability. Since we did not observe signatures of 
the Eckhaus reorganization in any of our simulations, we 
conclude that this mode of reorganization is not a generic 
mechanism in the Elastic Network (EN) model. 



ZZ reorganization of OD stripes during continuous 
area increase 

In our manuscript, we analyzed growth-induced ZZ 
reorganization of OD stripes in the EN model for in- 
stantaneous area increases. Here we point out that the 
characteristic sinusoidal bending of OD domains during 
this universal mode of reorganization is also observed for 
stripe-like OD pa tterns subject to continuous area in- 
crease ( Fig. 13). We initialized our simulations with 

sin(/cn 



a sinusoidal pattern o(x) ~ sin(/cmax^)- After a brief 
relaxation period of lOr we linearly increased the lin- 
ear extent of the simulated regions by 20%, 40%, 60% 
and 80% between t = lOr and t = llOr (correspond- 
ing to isotropic total area increases of 44%, 96%, 156% 
and 224%). Columnar layouts generally showed reorga- 
nization on similar timescales and the reorganization was 
more irregular when compared to the scenario of instan- 
taneous area increases (Fig. |13| :;, d). In all simulations, 
column spacing A exhibited a transient increase before 
the onset of reorganization followed by a decreased to- 
wards the initial value (data not shown). Handedness a 
decreased after onset of reorganization (data not shown) 
as expected from the characteristic sinusoidal bending 
of OD domains (Fig. [l3|. Note that stripe-like OD 



patterns subject to continuous area increase provide an 
intermediate scenario between OD stripes subject to in- 
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Adding columns at the periphery 
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Figure 11: Adding columns at the periphery requires large-scale remodelling. Schematic cross-sections through the cortical 
surface at an early postnatal stage (upper row) and an adult stage (lower row) are sketched. Dark/ white circles denote regions 
with ipsilateral/contralateral eye preference. Note that the number of neurons and the spacing of OD columns is assumed to 
be constant. In regions labelled with a red number, OD is required to change in order to preserve the initial OD layout in the 
central part. Note that this scenario implies large-scale rearrangements, in particular at peripheral regions of VI, and for some 
neurons multiple changes in OD over time. 
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Figure 12: Eckhaus reorganization upon large isotropic area increase (a, b) or area decrease (c, d). Upon area increase/decrease, 
OD stripes recover their initial spacing Ai by insertion/removal of OD stripes, (a) During isotropic area increase, the OD 
pattern transiently acquires a spacing A2 = Ai + ^A, with ^A > 0. By insertion of an new stripe the system recovers a spacing 
A3, roughly equal to Ai. (b) In Fourier-space, the initial pattern is centered at ±/cmaxX. ^A > corresponds to a shift of 
6k < 0. The initial spacing is recovered by the growth of new Fourier modes close to modes of the original pattern ±/cmaxX. 
(c) During isotropic area decrease, the OD pattern transiently acquires a spacing A2 = Ai + ^A, with 6A < 0. By removal of 
stripe the system recovers a spacing A3, roughly equal to Ai. (d) In Fourier-space, ^A < corresponds to a shift of 6k > 0. 
, with ^A > 0. The initial spacing is recovered by the growth of new Fourier modes at the center of the original pattern, i.e. 
i/cmaxX. Note that insertion/removal of columns requires strong area increase/ decrease, a regime not expected to be relevant 
during continuous growth. 



stantaneous area increase and realistic growth scenarios 
for complex OD patterns similar to the OD layouts in 
young kittens. 



Different Scenarios for Realistic Cortical Growth 

We tested three different simulation protocols toap- 
proximate the postnatal growth of cat VI (see Fig. |14^ , 
b): (i) linear increase of the linear extent with slope a, 
starting from Lq, 



L{t) = Lo^at, 
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Figure 13: ZZ reorganization of a stripe-like OD pattern in EN model during continuous isotropic area increase, (a, b) ZZ 
reorganization far from Turing Instability (r = 0.18, r] = 0.025, total area increase 196% between t = lOr and t = llOr). Until 
t = 50r the OD stripes remain almost unchanged. Subsequently, ZZ reorganization is apparent from the sinusoidal bending of 
domains (a) and the growth of Fourier Modes at /cmaxX ± qy (b) (c) ZZ reorganization close to Turing Instability (r = 0.04, 
other parameters as in a). Soon after the onset of ZZ reorganization the pattern acquires an irregular layout, (d) Power 
spectra for the OD layouts in c. Note, that despite rather complex reorganization, weak signatures of ZZ modes can still be 
identified. 



(ii) sublinear increase of the linear extent L, such that 
the simulated area A increases linearly with slope /3, 

L{t) = ^Ll + (3t, 

(iii) and logistic increase [56j of the linear extent L 

In every protocol, we increased the linear extent of the 
simulated region by a fixed amount (20%, 40%, 60% and 
80%) between t = Or and t = lOOr This corresponds 
to isotropic total area increases of 44%, 96%, 156% and 
224%. Note that according to [3j the area of cat VI in- 
creases between week 1 and week 12 by approximately 
150% percent. Simulations of OD development with dif- 
ferent growth protocols yielded qualitatively similar re- 
sults. Linear and sublinear increases described virtually 
identical growth scenarios ( Fig. |14^ , b) and provided 
a reasonable fit to the data on VI growth shown in [3j. 
Our systematic quantitative characterization of OD de- 
velopment in the EN model with growth (Figure 7, 8, 9 



in the manuscript) was carried out based on the scenario 
of linear increase of the linear extent (red traces in Fig. 
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Figure 14: Three different simulation protocols of realistic cortical growth, (a-b) Time courses of the linear extent L/Lq 
(a) and the area sizeA/Ao (b) for linear increase of area (green traces), linear (red traces) increase of the linear extent and 
logistic growth (blue traces). Simulations of OD development for all three growth protocols yielded qualitatively similar time 
courses of column spacing A, number of hypercolumns Nhc and bandedness a. Linear area increase and linear increase of the 
linear extent are virtually identical and both provide a reasonable fit to the data on cat VI growth in [3J. Our quantitative 
characterization of OD development in the EN model with growth (Figure 7, 8, 9) was carried out using linear increase of 
the linear extent (red traces). 



Videos 

Video 1. Development of OD column layout in the EN 
model without growth, starting from initial condition 
o(x,t = 0) = (77 = 0.025, r = 0.16). The video dis- 
plays OD maps (upper left) as well as time courses of the 
parameters mean column spacing A (upper right), num- 
ber of hypercolumns Nhc (lower left), mean bandedness 
a (lower right), characterizing OD layouts during devel- 
opment. At t = 4t OD columns have segregated. After 
lOr, columns start reorganizing and merge progressively 
towards a stationary stripe-like layout. Whereas A, and 
Nhc reach near mature levels around lOr the banded- 
ness a increases during the entire time course, capturing 
the slow reorganization of OD columns. 
Video 2. Zigzag reorganization in the EN model after 
instantaneous isotropic area increase. The video displays 
OD maps (upper left) as well as time courses of mean col- 
umn spacing A (upper right), number of hypercolumns 
Nhc (lower left), mean bandedness a (lower right). Af- 
ter instantaneous area increase at t = lOr, OD stripes 
bend sinusoidally, a behavior typical for ZZ reorganiza- 
tion. A jump in A is observed shortly after t = lOr. 
Subsequently, however, A decreases to almost the initial 
value and consequently, Nhc increases, a decreases, re- 
flecting the bending of OD domains. When increasing 
the system size, the EN model displays all features typi- 
cal for ZZ reorganization. 

Video 3. OD development in the EN model with growth 
at near realistic growth rate (156% total area increase 
between t = and t = lOOr, r] = 0.025, r = 0.16). 
The video displays OD maps (upper left) as well as time 
courses of mean column spacing A (upper right), num- 
ber of hypercolumns Nhc (lower left), mean bandedness 
a (lower right). Despite the large area increase, column 
spacing A increases only mildly and transiently, imply- 
ing a strong increase in Nhc- Note, that growth-induced 
reorganization of OD layouts shows aspects of the ZZ re- 
organization that is captured by a drop in a. 



